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We give a concise formula for the non-Gaussianity of the primordial curvature perturbation gen- 
erated on super-horizon scales in multi-scalar inflation model without assuming slow-roll conditions. 
This is an extension of our previous work. Using this formula, we study the generation of non- 
Gaussianity for the double inflation models in which the slow-roll conditions are temporarily violated 
after horizon exit, and we show that the non-linear parameter fML for such models is suppressed by 
the slow-roll parameters evaluated at the time of horizon exit. 
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I. INTRODUCTION 

Non-Gaussianity of the primordial curvature perturba- 
tion is a potentially useful discriminator of the many ex- 
isting inflation models jll,l2]- Planck [3| is expected to de- 
tect the primordial non-Gaussianity if the so-called non- 
linear parameter, Jnl, which parameterizes the magni- 
tude of the bispectrum, is larger than 3 ~ 5 [l|,|j|. Higher 
order correlation functions such as trispectrum would 
also be a useful probe of the primordial non-Gaussianity 
0) H 0| • Hence it is important to theoretically under- 
stand the generation of non-Gaussianity. 

Standard single slow-roll inflation model predicts 
rather small level of the non- linear parameter, JnLi 
suppressed by the slow- roll parameters @]. In multi- 
scalar field inflation models with separable potential, ex- 
plicit calculations in the slow-roll approximation show 
that /atl is also of the order of the slow-roll parame- 
ters [1 0, [m, m, in, [ii, [H, [ii . This feature is likely to 
hold for more general potential that satisfies the slow-roll 
conditions [13,[l3j[3- On the other hand, the generation 
of non-Gaussianity due to the violation of the slow-roll 
conditions remains an open problem. 

In this paper, based on 5N formalism [20, [2l|, [23, [23, 
[24| , we give a useful formula for calculating the non-linear 
parameter in the multi-scalar inflation models includ- 
ing the models in which the slow-roll approximation is 
(temporarily) violated after the cosmological scales exit 
the horizon scale during inflation. Current observations 
do not exclude such models. We also study the possi- 
bility of the large non-Gaussianity for double inflation 
model [2^, [2y| , as an example of such a model. 

In section [nl we derive the formula for the non- linear 
parameter Jnl, which is an extension of our previous 
paper [l9|. In section HITl we analyze the primordial non- 
Gaussianity in the double inflation model using our for- 
mula. We give a summary in section [TVl 



II. FORMULATION 



In this section, we derive a formula which is an exten- 
sion of that given in our previous paper [19| to the non 
slow- roll case. We use the unit Afpj — {%-kG)~ = 1. 



A. Background equations 

We consider a A/'-component scalar field whose action 
is given by 



S^ 



\hijg^^''d^cj,'d,^' + V{(^) 



, (1) 



(/,J=1,2,...,AA) , 



where g^^ is the spacetime metric and hij is the metric 
on the scalar field space. In the main text, we restrict 
our discussion to the flat field space metric hjj = 6jj to 
avoid inessential complexities due to non-flat field space 
metric. Extension to the general field space metric is 
given in appendix [Bj 

We define ipl{i — 1,2) as 



^(^ 



fi = 



dN 



(2) 



where dN = Hdt with H and t being the Hubble parame- 
ter and cosmological time, respectively. Namely, we take 
e-folding number, N , as a time coordinate. For brevity, 
hereinafter, we use Latin indices at the beginning of Latin 
alphabet, a, h or c, instead of the double indices, /, i, 
i.e.,X- = X/. 

Then, the background equation of motion for (p" is 



ip" = F^iifi) 



(3) 
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where F°'{— F.l) is given by 



Fi 



'fii 



Fi 



V_ 



I y 



(4) 



where V^ = S^'^{dV/d4>'^), and the Friedmann equation 
is 



H' = 



2V 



&-fif2i 



(5) 



with if2i = 5ijtf2 



(2) 



Evolution equation for (5<P° is given by 



d „(2), 



dN 



Sip^^N) = P'l{N)Sip''{N) 



+Q\,{N)S^\N)S^'{N) , (12) 



where Q^^c i^ defined by 



B. Perturbations 

In the SN formahsm [H, [H [H, the evolution of the 
difference between two adjacent background solutions de- 
termines that of the primordial curvature perturbation 
on super-horizon scales. In this paper, we use the word 
"perturbation" to denote the difference between two ad- 
jacent background solutions. In this subsection, we an- 
alyze the time evolution of the perturbation and relate 
the result to the curvature perturbation. 

The solution of the background equation ^ is labelled 
by 2A/' integral constants A"^. Let us define dip"" as the 
perturbation. 



Sip^iN) = ip''{X + 6X; N) - </?°(A; N) 



(6) 



where A is abbreviation of A" and SX"' is a small quantity 
of 0(15). 6(p°'{N) defined by Eq. ([6]) represents a pertur- 
bation of the scalar field on the N — constant gauge [23] . 
For the purpose of calculating the leading bispectrum 
of the curvature perturbation, it is enough to know the 
evolution of S(p°'{N) up to second order in 5. For later 
convenience, we decompose Sip"^ as 



(1) 1 (2) 



(7) 



(1) (2) 

where Sv"" and Sv"" are first and second order quantities 
in S, respectively. 



Evolution equation for Jv" is given by 



dN 



Sv^iN) = P%{N)Sv''{N) 



ij\ 



where -P'^(= PiJ) is defined by 






(0) 



(8) 



(9) 



(0), 



Here <^(7V) represents the unperturbed trajectory. The 
explicit form of P^j, is shown in appendix [X] Formally, 
solution of this equation can be written as 



5^%N) ^ K\{N,N,)6v\N,) , 



(10) 



where K\ is a solution of 



Q\ 



d^pa 



be 



dip^dip^ 



dPl 



(0) 



dif"^ 



(0) 



(13) 



The explicit form of Q\ci~ Quk) ^^ shown in appendix 
\X[ Let us choose the integral constants A" as the initial 
values of ip"" at N — N^,, namely, A° — ip'^{N^,). Then we 

have Sip"{N^) = (5A". Hence S^^iN) vanishes at N^. Un- 
der this initial condition, the formal solution of Eq. (fT2|) 
is given by 



5<^'°(7V) 



N 



N, 



dN'A^,{N,N')Q'^^{N') 

x6'^^''{N')6^^^'^{N') . (14) 



According to the 5N formalism [22, llj], the curva- 
ture perturbation on large scales evaluated at a final 
time, N = Nc, is given by the perturbation of the e- 
folding number between an initial flat hypersurface at 
N — Nf, and a final uniform energy density hypersurface 
at A^ = A^c- Let us take A^* to be a certain time soon after 
the relevant length scale crossed the horizon scale, H~^, 
during the scalar dominant phase and Nc to be a certain 
time after the complete convergence of the background 
trajectories has occurred. At A^ > A^c the dynamics of 
the universe is characterized by a single parameter and 
only the adiabatic perturbations remain^. Then, the e- 
folding number between TV* and Nc can be regarded as 
the function of the final time Nc and (p°'{N^,), which we 
denote NiNcipiN^)). 

Based on SN formalism, the curvature perturbation 
on the uniform energy density hypersurface evaluated at 
N — Nc is given by 



C(A^e) ^ 5N{Nc,ipiN,)) 

Jipt + IrNabJip'^Sip'i 



N,, 



(15) 



where 5(p'^ = d(p'^{N^) represents the field perturbations 
and their time derivative on the initial flat hypersurface 
at A^ == A^,. The left hand side in Eq. ^ is obvi- 
ously independent of the initial time A^*, and hence so 

is 6N{Nc,ip{N^)). Here we also defined Na* = NaiN^,) 



—A%{N,N') ^ P'^c{N)A%{N,N') , 
with the condition A%{N, N) = Afj(Ar, A^) = S^jS^ 



(11) 



^ We consider the case in which isocurvature perturbations do not 
persist until later. 



and Nab* = Nab{N^) by 



NaiN) 
NabiN) 



_ dN{N,,ip) 



dip'' 



(0) 



dip'^'-dtp'' 



(0) 



(16) 

(17) 



evaluated at N = N^,. 

It is well known that the curvature perturbations on 
an uniform density hypersurface, (^, remain constant in 
time for N > Nc- Hence, C(-^c) gives the final spectrum 
of the primordial perturbation. 

Let us take Np to be a certain late time during the 
scalar dominant phase. Then we have 



({Nc 



SN{N,,^{Nf)) 
1 
2 



NapS^^F + ^NabpS^^pS^], + ■ ■ ■ . (18) 



where 6ip}, = (5(^'^(7Vf), NaP = NaiNp) and NabP = 
NabiNp). During the period with N^, < N < Np, we 
can use the solutions for 6ip'^ given by Eqs. pU|) and pi)) . 
Using these solutions, we obtain the relations; 



Na,^NbFA\iNF,N,) , (19) 

Nab* = N,dpA\{NF,N,)A%iNF,N,) 

i-Nf 

+2 / dN'N,{N')Q'a^iN') 

xA'^^{N',N^)A%{N',N,) , (20) 

Na{N) = NbFA\{NF,N). (21) 

C. Non-linear parameter 



with 



In this subsection, we derive a formula for the non- 
linear parameter /jvl by making use of the SN formal- 
ism. We first give the definition of /nl- It is defined as 
the magnitude of the bispectrum of the curvature per- 
turbation c. 



Bi^(ki,k2,k3 



6 fwL 

5 (27r)^/^ 



Pdki)Pdk2) 



+Pc{k2)Pdk3) + Pc{k3)Pc{ki) 



(22) 



where P(^ is the power spectrum of (.- The definitions of 
P^ and B(^ are, respectively, 

(Ck,Ck.)=<5(ki+k2)Pc(fci) , (23) 

(CkiCk.Ck3> ^ <5(ki +k2 +k3)Bc(fci,fe,fc3) . (24) 



Equation ([22]) restricts the form of the bispectrum. The 
bispectrum in general does not take that simple form. 



In fact, sub-horizon perturbations of fields give differ- 
ent fc-dependent form of the bispectrum [17|. However, 
the sub-horizon contribution to the bispectrum is sup- 
pressed by the slow-roll parameters evaluated at the time 
of horizon exit^. In contrast, the super- horizon evolution 
always gives the bispectrum in the form of Eq. (j22p inde- 
pendent of the number of fields (see below). If /nl ^ 1, 
which is an interesting case from the observational point 
of view, then the contribution due to super-horizon evo- 
lution dominates the total bispectrum. 

We assume that the slow-roll conditions are satisfied at 
N ~ N^:. Then, to a good approximation, <5(p{^ becomes 
a Gaussian variable [8|, ll7J | with its variance given by 



is^ij^t) = 6" ^^y 



(25) 



and (^2 becomes function of (p^. Differentiating (^2 



— -^, we have 



Sfi 



2* 



V'Vj V' 



y2 



V 



s^(* 



(26) 



The higher order terms are also suppressed by the slow- 
roll parameters. Hence, ^1^2* i^ Gaussian as is 6ipl^ to 
a good approximation. Then, we can write down the 
variance of Sip"^ as 



{S^iSifil 



A 



ab 



2^ 



(27) 



At the first order both in the field perturbation and slow- 
roll limit, the matrix A'^'' 



A\f can be written as 



aU _ /r/J aU — aU — f 

^11 ~ " I ^12 ~ ^21 ~ ' 



IJ 



A 



IJ 



■ K'^ 



KJ 



where 



JJ 



y^(0)y-^(0) Y^-d<t>) 



Since e 



/,/ 






Vict,) 



(0) 



(28) 



(29) 



: find that {6ip{^5Lp2*) and 
{5f2*^f2*) ^'^s suppressed by the slow-roll parameters. 
At the same order, 0{e,r]), it is known that we need to 
add the slow-roll correction terms to A{i in Eq. 
Such corrections to A{f are given in |28|]. 



^ Here, we consider the case that inflation is induced by the canon- 
ical scalar fields. In such a case, current observations of the spec- 
trum of curvature perturbation restrict the models so that the 
slow-roll conditions are satisfied until the cosmologically relevant 
scales exit of the horizon scale. Otherwise, for example, in Dirac- 
Born-Infeld (DBI) inflation model, the sub-horizon contribution 
to the bispectrum not only has different fc-dependent form, but 
also can be large enough to be detectable in the future experi- 
ment without inconsistency with current observations (27ll . 



Using these equations, to the leading order, the non- 
hnear parameter is written as 



1 



{Na.&i) 



NabFe''{NF)Q\Np 



Nf 



+ / dN'N,iN')Q^,,iN') 



xe''iN')e''{N') 



where 



Q^iN) = A^^iN, N,)A''''Nb. 



(30) 



(31) 



and 8" = 0"(A^*). As we mentioned before, we have 
neglected the non-Gaussianity from the sub-horizon con- 
tributions in deriving Eq. (|5(I)l . Eq. ([5D|) shows that, aside 
from NaF and NatF whose explicit form will be given in 
the next subsection IIIDl Jnl is completely determined 
by the quantities Na{N) and <d"'{N). These quantities 
obey the following closed differential equations. 



dN 
d 

dN 



NaiN) = ~NbiN)P''^{N) 
e^iN) = P'l{N)e''{N) . 



(32) 
(33) 



First, we solve Eq. (|32|) backward till N ^ N^, under 
the initial conditions NaiNp) = NaF- Then we solve 
Eq. ((33|) forward till N = Np under the initial condi- 
tions Q°-{N^) = A°■''Nl,.^.. Substituting these solutions 
into Eq. (PU)) . we obtain Jnl- 

The evolution equation of 6° (TV), Eq. ^, is iden- 
tical to that of Sip"", Eq. ([8]), at the linear level. Since 
Eqs. ([5^ and (P5|) are mutually dual, a variable com- 
posed of Na{N^) and e"(iV,); 



W{N,) = Na{N,)Q''{N,) , 



(34) 



becomes constant irrespective of N^. This constancy of 
W{N^:) corresponds to the constancy of SN{Nc,(p{N^,)) 
which was mentioned below Eq. (jl5[) . 



Here, let us evaluate C{Np) on the uniform energy den- 
sity hypersurface at iV = Np, neglecting the later evo- 
lution of the curvature perturbations during the period 
with Nc > N > Np. In this case we can obtain explicit 
forms of NaF and NabF, which appeared in Eqs. ((2T|) 
and ()30|) . written in terms of the background quantities 
at A^ = Np. On the super- horizon scales, the uniform 
energy density hypersurface is equivalent to the constant 
Hubble hypersurface. Then, C(-^f) is evaluated by the 
time shift SN measured from the H = constatnt hyper- 
surface. Therefore at A = Np we have the equation; 



H [ip'^iNp + C(Af))) = H{v>''{Nf)) 

The Hubble parameter H is given by Eq. ([ 
Eq. ([55]) with respect to ({Np), we obtain 



with 



CiNF) 



Nal 



NapSifl -f -NabF^iflSip^p, 



Haiv) 



NabF = - 



Uabi'P) 



(0) 



H,iip)F-{ip) 



(0) 



(35) 
Solving 

(36) 

(37) 
(38) 



where 



Uab — Hab — j-j- rn^ (HcP'^b 



Hdpd 

2HaHb 



F^HrAF" 



F'Hab) 

+ F'PjHa) 



(39) 



and Ha = dH/dip", Hab = d^H/dip'^dip''. The explicit 
forms of Ha and Hab are shown in appendix [Xj The 
right-hand sides of Eqs. ([57)1 and ([55]) are written in terms 
of local quantities, i.e., those depending only on 0^ and 
d(p^ /dN a,t N — Np. Hence, once we specify a cen- 

(0) 

tral background trajectory </5°, we can readily determine 
NaF, NabF- SN evaluated using the determined NaP 
and NabF depends on Np unless Np > Ac, where Nc is a 
time when the complete convergence of the background 
trajectories occurs. 



D. Non-linearity generated in scalar dominant 
phase 

In order to evaluate Eq. ([50)) . we need to know NaP 
and NabF - If one takes into account the evolution of the 
curvature perturbations after the scalar dominant phase, 
isocurvature perturbations may remain during preheat- 
ing/reheating era after inflation. In such cases, in order 
to calculate NaP, NabF, we need to investigate the evolu- 
tion of the background e- folding number with the effect of 
short wavelength/radiation component, which is beyond 
the scope of the present paper. 



III. DOUBLE INFLATION MODEL WITH 
LARGE MASS RATIO 

Substituting Eqs. §7^ and dSS]) to Eq. jSO]), the non- 
linear parameter, /nl, can be also defined as a function 
of the final time, Np. As mentioned above, the curvature 
perturbations on uniform energy density hypersurface, C., 
evolve when Np < Nc, and hence /atl also evolves. 

In this section, using the formulation given in the pre- 
vious section |TT1 we calculate the non-linear parameter 
fNL{Np) for the double inflation models, which violate 



the slow-roll conditions for a certain period when the cos- 
mologically relevant scales are well outside the horizon. 
The potential is given by [23, [2^ , 



Vi(t>,x) = 2™^*^^ + 2™^'^^ 



(40) 



We assume a large mass ratio, i.e., m^/m^ » 1. Because 
of m^/rricf, ^ 1, the energy density of x-field, p-^, decays 
faster than that of ^-field, p^. Hence if p^ ^ p^ at initial 
time, py^ never dominates the energy density during the 
later evolution of the universe. In such a model, only a 
single chaotic inflation induced by c/i-field occurs. It is 
known that in a single scalar slow-roll inflation, /jvl is 
suppressed by the slow-roll parameters. Here we assume 
the opposite, p^ S> p^, at the initial time. In this case 
the slow-roll conditions are badly violated when (f> — x 
equality (p^ ~ p^) occurs. We assume that both fields 
are in the slow-roll phase at the initial time and the in- 
flation induced by 0- field also occurs after — x equality. 
Denoting the initial values of fields at N — N^, as 0* and 
X*, this assumption implies (f>*,x* ^ 1- 

There are several works in which the primordial non- 
Gaussianity in the two-scalar chaotic inflation model 
whose potential is given by Eq. (|1(I)) was investigated. 
In Ref. [9|, Il0l | , the authors analyzed the primordial non- 
Gaussianity generated during A/'-flation 29] model and 
gave a simple analytic formula using the slow-roll ap- 
proximation even on super-horizon scales based on SN 
formalism. In Ref. [i3| , using their slow-roll formula the 
authors also calculated the non-linear parameter /jvl for 
the two-scalar chaotic inflation model whose potential is 
given by Eq. (j40p with the mass ratio m^/m^ = 9. In 
such case, the slow-roll conditions are not violated during 
scalar dominant phase, so they could use the slow-roll 
formula. They compared the result obtained by using 
their analytic formula with that obtained numerically, 
and they found that the non-linear parameter does not 
become large in such case. In Ref. ^i3], the authors also 
provided another formulation for the non-linearity gener- 
ated on super-horizon scales without slow-roll conditions. 
They also calculated the non-linear parameter using their 
formulation for the potential given by Eq. (j40p with the 
mass ratio m^/m-)^ = 12. They also found that the non- 
linear parameter evaluated at the end of inflation does 
not become large in such model. 

Here, using the formulation given in the previous sec- 
tion, we numerically calculated /nl for the mass param- 
eters 7710 = 0.05, my^ 1.0 and compared the result with 
the previous work |16J. 



A. Numerical calculation for m^/m^ = 20 

Using the formulation given in the previous section, 
we numerically calculated /nl for the mass parameters 
m^ = 0.05, m^ = 1.0. Hence the mass ratio is m^/m^ ~ 
20. Initial value of flelds are set to 0* = x, = 10. We 
choose the e-folding number as a time coordinate and 



set the initial time TV* to 0. The central background 
trajectory in the field space is shown in Fig [T] and the 
evolution of these two fields as a function of N is shown 
in Fig [21 We define slow-roll parameters as 



i]<i><i> 



1 dH 
H'dN 

'' V ' 



%x 



V 



(41) 
(42) 



The evolution of the slow-roll parameters is shown in 
Figs. [1[1 and [H 








FIG. 1: This figure shows a background trajectory in field 
space. We set initial values of fields to (j>f = x* ~ 10. 

From Fig. [21 we see that initially x-^eld decays rapidly 
while 0-field remains almost at its initial value due to 
large Hubble friction. At around TV = 26, x-field reaches 
zero and starts damped oscillation. As shown in Fig. [3l 
around this time, the slow-roll parameter e exceeds 1 for 
a moment and the infiation ends once. From Figs. [Hand 
[5l also 77^^ exceeds 1 after x-field settles to the minimum, 
while, (/)-field is slow-rolling during this phase. 

In calculating the non- linear parameter /atl, we regard 
Jnl as a function of Np in the same sense as C,{Nf) in 
Eq. pop . We show the evolution of the non-linear param- 
eter fNLiNp) in Fig- [3 From this figure, we find that 
the non- linear parameter fNLiNp) temporarily oscillates 
with the maximum amplitude reaching ~ 8. As soon as 
X-field settles to its minimum, the universe becomes dom- 
inated by ^-field due to the decay of p^ , and the universe 
starts the second inflation driven by (^-field. After the 
second inflation starts, iso-curvature perturbation decays 
rapidly due to the decay of x- Hence the time Nc is before 
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FIG. 2: (color online) This figure shows evolution of (j> (dashed 
red line) and x (solid blue line). We choose the e- folding 
number, A'', as a time coordinate and set the initial time, A^*, 
to 0. 



FIG. 5: This figure shows the evolution of the slow-roll pa- 
rameter r)^-^ defined by Eq. (|42p . 
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FIG. 3: This figure shows the evolution of the slow-roll pa- 
rameter t defined by Eq. (|4ip . 
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FIG. 4: This figure shows the evolution of the slow-roll pa- 
rameter r]^tt> defined by Eq. (|42|) . 
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FIG. 6: This figure shows the evolution of non-linear param- 
eter, fNLiNp)- After the second infiation, Nf > 30, the 
background trajectories in phase space have completely con- 
verged. Thus, the curvature perturbation on a constant Hub- 
ble hypersurface, ^(Nf), remains constant on large scales at 
Nf > 30 and the /jvi(A^F) also remains constant in this era. 



the end of the second inflation in this model. Then, the 
sum of fMhi^F > 30) and the sub- horizon contributions 
gives the primordial non-Gaussianity, independent of the 
reheating process. The final value of fNLiNp > 30) is 
0.01004. Hence the generation of large non-Gaussianity 
due to the violation of the slow-roll conditions in this 
model does not occur. 

In Ref. [ig|, the authors numerically calculated the 
non-linear parameter for the model whose potential is 
given by Eq. (|40p with mass ratio m^/m^ — 12. They 
used their own formalism, which is different from SN 
formalism used in this paper. The behavior of the non- 
linear parameter /jvl obtained by their calculation seems 
similar to our result shown in Fig. [6l 



B. Approximate analytical expression for /a 



In the previous subsection (jIII All , by numerical cal- 
culation, we found that the final value of /nl is much 
less than 1, although the curvature perturbation becomes 
highly non-Gaussian {Jnl — 8) for a moment. We can 
derive the analytical expression for the final value of /nl 
in the limit of large mass ratio m-^/m^ ^ 1, which is the 
subject of this subsection. 

The background equations are given by 






r2 d(j> , _2 
dN 



3H'^+mlcb^Q, (43) 



3i?^^+m^X = 0, (44) 



H 



2 _ ml(p^ + m\^ 



dN 



dN 



(45) 



Within the slow-roll approximation, these equations re- 
duce to 



Here we recall that we have set A^* = 0. After the time 
when X ~ 1, ^-field is no longer in the slow- roll phase 
and p^ evolves as; 

p^{N) ~ p^{x{N,,))e^^[-i{N ^ N,,)] , (49) 

where Nos represents the time when ^-field starts oscil- 
lation. From Eq. ([35]), we have 



A^„ 



4 



x{No 



Then, Eq. 



can be rewritten as 



Px(^)-Pxexp(-3A^+-X: 



(50) 



(51) 



where p^ is independent of A^, x* and 0*. Since p^ de- 
cays in time for A^ > Nos , P<p dominates the total energy 
density of the universe at some time and (^-field starts 
slow-rolling. Let us denote the e-folding number at the 
time when p^ ~ p-^ as A^eg- For A^ > Nf,q, the evolution 
of ^-field is given by 

<I>{N) ~ ^q^l-4{N-Ne,) , (52) 



A^eq is obtained from the 



during the slow-roll phase, 
equation, 

P^iNeg) ^ -ml<t>{N,qf ~ p^{N,q) . 
Combining Eq. (|?T|) with Eq. ([55]) . we have 



Y 1 



2Px 



(53) 



(54) 



The total energy density after A^ — Neq is approximately 
P0. Hence, using Eqs. ([5^ and ([5^ . we have 



d 

Za' 



6 






dN' 

2 2\ 

m^X ) 



mix 



(46) 



At the initial time, N = N^,, both fields are in the slow- 
roll phase and </)* ~ %* » 1 . Then the ratio of the time 
derivatives of the scalar fields is 



piN) ~ p^iN) 



'-^\N) 



0, + X* - 4 A^ + 



1 , ™'<^2 



(55) 



d(l)/dN ^rn^ff) 
dx/dN m?, X 



<1 



(47) 



and also p^ 3> Pc/, in the limit of large mass ratio. Hence, 
to a good approximation, 0(A^) ~ 4>* until the energy 
density of 0-field dominates the total energy density of 
the universe, i.e., p^ ~ p^- On the other hand, solv- 
ing Eq. ([46]) , the evolution of x- field during the slow-roll 
phase can be obtained as 



X(A^) ~ Vxl - ^N 



(48) 



After A^ — Neq, the universe is dominated by the sin- 
gle slow-roll component (/)-field, and the curvature per- 
turbation on the uniform density hypersurface becomes 
constant in time on super- horizon scales. From Eq. (j55p . 
we obtain an expression for the e-folding number at the 
final time in terms of </>* and x* as 



N{Nf,(J),,X*) 



XZ) - olog' 



PF 



1, ^0 
- log — — 
3 ^2px 



(56) 
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where pp = p{Np). This gives 



5 



NjNjN" 
2 



4>i + xi 



i + o{ 



')) , (57) 



where TV/ = dN/dcjji, Njj = d^ N / d(l)id(j)i ( Here, c/)^ = 
((/), x)). Hence Jnl is 0((/)^^, X^^)i and it is much less 
than 1. Jnl estimated using this expression is w 0.01 for 
m^ = 0.05, m^ = 1.0 and (j)* = X* = lOj which agrees 
well with the numerical result in the previous subsection 
(seeHIEl. 

The reason why we get small /nl in the double infla- 
tion model can be understood as follows. In the case of a 
large mass ratio, evolution of the universe can be clearly 
divided into three stages. First one is the inflation in- 
duced almost by a single field, i.e., x-field. The second 
one is the non-inflationary phase where the energy den- 
sity of the oscillating x-field is still larger than that of the 
(^-field. The third one is the inflation induced by a single 
field, i.e., 0-field. Then the total e-folding number during 
all these stages mainly comes from the two inflationary 
stages, i.e., the first and the third stages. The perturba- 
tion of the e-folding number generated during the second 
stage is negligible. Hence as a crude approximation, the 
dynamics can be regarded as just the sum of two single 
field infiationary phases. This picture gives the first term 
of Eq. (j56p . The correction to this picture is represented 
by the second term, which is indeed minor for A^ ~ 60. 

The derivation of Eq. ([57)1 can be straightforwardly 
extended to more general double inflation models with 



y(^,X)-Ci02P + c2X 



2p 



(58) 



where p is an arbitrary positive integer. In this case, the 
final value of Jnl is more suppressed than Eq. ([57|) by a 
factor 1/p. 



C. Comment on /jvl in A/^-flation model 

Discussion of the previous subsection about the small 
/jvi in the double infiation model gives us some insight 
into the generation of the non-Gaussianity in the so- 
called A/'-flation model [29], in which the potential is 
given by 



n0) = Ef^^^- 
1=1 



(59) 



This is a generalization of the double inflation model to 
an arbitrary number of fields. We consider the case in 
which all fields have large initial amplitude, i.e., (/<^ 3> 1. 
We also assume m^ < m^ < ■ ■ ■ < m^ . In this case, 
(fr field decays first. If pj^ dominates the total energy 



density of the universe until (j)^ ~ 1 ^, then the inflation 
is almost like the single field inflation by </> -field. After 
(/)-'^-field decays, then 0-™"^ -field decays and so on. Hence 
the leading e-folding number would be given by 



TV: 



1 



7=1 



(60) 



Then the corresponding f^L is 



JNL 



<S^'T'~^- 



(61) 



Roughly speaking, /jvl is suppressed by the inverse of 
the e-folding number. This result is quite similar to that 
in Rcf. [g| , where the authors studied under the slow-roll 
approximation. 



IV. SUMMARY 

Based on the 5N formalism, we have derived a useful 
formula for calculating the primordial non-Gaussianity 
due to the super-horizon evolution of the curvature per- 
turbation in multi-scalar infiation without imposing slow- 
roll conditions. This formula can apply for the infla- 
tion models with general field space metric, hjj^ as far 
as super-horizon contributions are concerned. Generally, 
when one calculates the non-Gaussianity of the curvature 
perturbations, one has to solve the second order pertur- 
bation equations. In doing so for a multi-scalar infiation, 
there appear tensorial quantities with respect to the in- 
dices of the field components. Our formula reduces the 
problem of calculating the non-linear parameter /jvl to 
solving only first order perturbation equations for two 
vector quantities. This reduces 0{JV^) calculations to 
0{J\f) ones where M is the number of the scalar field 
components. Hence our formalism has a great advantage 
for the numerical evaluation of Jj^il in the infiation model 
composed of a large number of fields. 

We have also studied the primordial non-Gaussianity 
in double inflation model as an example that violates 
slow-roll conditions by using our formalism. We found 
that, although f^L defined for the curvature perturba- 
tion on a constant Hubble hypersurface exceeds 1 for a 
moment around the time when the slow-roll conditions 
are violated, the final value of Jnl is suppressed by the 
slow- roll parameters evaluated at the time of horizon exit. 
We have shown that this can be understood even analyt- 
ically in the 5N formalism. This result is straightfor- 



^ If PaA becomes subdominant before t/i-™ ~ 1, then <j/^ is slowly 
rolling until H decreases to mj^. Such a case was studied in 
13 llOll under the slow-roll and horizon crossing approximation. 
The final value of Jml was found to be suppressed by the slow- 
roll parameters. 



wardly extended to more general double inflation model 
and A/'-flation model. 
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APPENDIX B: EXTENSION TO THE GENERAL 
FIELD SPACE METRIC 

In this appendix, we will extend the formulation given 
in sec. II for the flat field space metric hjj = Sjj to the 
general case. For general case, the background equations 
corresponding to Eq. ^ for flat field space are given by 



dN ^ ' dN 



^2 ' 



(Bl) 



APPENDIX A: SPECIFIC EXPRESSION 

In this appendix, we show the explicit forms of P\{^ 
P/i) in Eq. ®, QV(= QIjk) in Eq. (USD, i/.(- H\) in 
Eq. 1211) and Hab{= H'lj) in Eq. ^. P-j is given by 



p/i 
pii 

^2. J 



pI2 



s', 



V 



V 






P2.J = V2V2J + -y'P2j - -^5 J 



(Al) 



where DA^ = dA^ + T^ jj^A-^ dip^ , T^jj^ is Christoffel 
symbol with respect to the field space metric hu, and 



Fi 



vi 



Fi 



V 
IP 



V2 



h 



V 



(B2) 



In the same way as Eq. 
tion as 



we can define the perturba- 
(B3) 



S^^iN) = (^"(A + 5\- N) - ip^iX; N) 
and we decompose S^pf as 



(1) 1 (2) 

Sip"" = Sip" + -Sf" 



(B4) 



Q^^by 



The curvature perturbation can be expanded in terms of 
S(pj as 



Qi 



111 



1/22 



JK — QlJK — QlJK — ^IJK 



Q 



^ = 



^2JK 




y'^yi 



jVK 



V v^ 

y'xyj y'yjK 



qI12 



V 
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2, IK 



V 

~v 



K 



1/2 

V'Vj 

y'yx 

y2 



V2 



V2K 



^21 



yyjyi 



K 



y3 



qI12 
W2KJ 



Q2JK = " K^2J + jP2K + OjK "• 



Y! 

V 



¥>^ + — . (A2) 



H^j is given by 



H] 



HVi 



Hf 



IP 
2V 



•^21 



(A3) 



C ~ Na5p'' + -Nab5v''5v^ + ^ 



(B5) 



However, Sipl defined in this way does not transform as 
a vector under the "coordinate transformation", (j)^ -^ 
(j)^ {(j)) because Eq. (jB3p takes the difference of variables 
at different points in field space. Hence the evolution 
equation for S(pl lacks the manifest covariance. To avoid 
this drawback (though not essential), we introduce new 

(1) 
perturbation variables Sp"" = S(p°' 



(2) 

^Sp"" defined by 



Sp>{ 



(1), 

S<P2 



— a, 5p,,-. 

^^2.. Pi 
^a, 5p,, 



D dp}{ 
dX~dX 

_ DW2 

' dA2 






(B6) 

(B7) 



which transform as vectors. Differentiating the back- 
ground equations (jBip with respect to A, we obtain the 
evolution equations for the perturbations of the scalar 



10 



fields in curved field space metric. To the first order, we 
have 

^4'"(A^) = iP%iN) + AP%{N)) 4\iV) , (B8) 



where the P°j, is given by 



DF" 



pa __ 



(0) 



(B9) 



and AP"^, which vanishes in the flat field space metric 
representing terms due to the curvature of the field space, 
is given by 



API) = AP/j = AP/j = 
AP/i = -P'l,7K^2 ^2"" ■ 



(BIO) 



Here, R^ j^jx represents Riemann tensor associated with 

(2) 

the field space metric hij . Evolution equation for 5(p'^ is 
also given by 

§;^Sr{N) = {P%{N) + AP%{N)) 6^\N) 

+ iQ%AN) + AQ%,{N)) d^'>{N)S$%N) , 

(Bll) 



where 



of the first and second order equations are given by 



S^''{N)^A%{N,N,)S^\N,), 

(2) fN 

50''{N) = / dN'A\{N,N') 

JN, 

X [Q\AN') + AQ\^{N')\ 4^{N')4\N') . 

(B15) 
The expression for the curvature perturbation in terms 
of 5tf^ , Eq. (|B5[) , can be also rewritten as that in terms of 
5(p°' . To make the expression more concise, we introduce 
new variables Na and Nab as expansion coefficients of C 
in terms of 5(^° so that we have 



C ~ NaF5^''F + ^NabF5^'i'S(p^F + ' 



(B16) 



Comparing the two expressions (|B5[) and (|B16|) . we find 

Nj = Nj - iV|.r5,^2^ , Nf = Nf , 
N}} = Nj} + Nf,T\,rT\,^i'^^ 

-N],ViTljcp^ - NJ,T\j + N],T^j^T\jcp^ , 



Njj = NPj - iv|.r^^ - ^|?jr^i^^ 

N]} = Nfj - NlY^.j - N]i,T^j^v'i 
Nfj = Nfj . 



(B17) 



D^F° 



be 



dif'^dip'^ 



DPI 



(0) 



dip'^ 



(0) 

ip=V(N) 



,(B12) 



and the explicit form of AQ\^ is given by 

AQijK = -R JKLV2 ) 

AQifx = AQifx = AQifx = , 

-2P KJL'P2 : AQ2JK = -R LK.jf2 : 



AQ2JK 

AQ 
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2JK 







(B13) 



where V k is covariant derivative with respect to hjj. 
By introducing the matrix A as the solution of 



From the constancy of NaSip"", we find 

Na{N) = NbFA\{NF, N) . (B18) 

We also redefine 9f in Eq. ([3T|) as 

e'^(7V) = A%{N,N^)A'"'N„ , (B19) 

where A"-'' is defined by 



~a T~b \ 



{S^tSif 



A 



ab 



2^ 



(B20) 



D 



—A%iN,N,) = {P'^,iN)+APUN))A^,iN,N,), 



A°'^ for general case also depends on the field space met- 
ric, hij. 

Then, the expression for the non-linear parameter in 
general cases is also given by Eq. ([5D|) with Nc{N') and 
with the condition A[j(A^, N) = 5'j5.\ formal solutions Q\bi^') replaced by Nc{N') and Q\^^{N') + AQ''^^{N'). 



(B14) 
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